Topological dynamics and dynamical scaling behavior of vortices in a two-dimensional 

XY model 
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By using topological current theory we study the inner topological structure of vortices a two- 
dimensional (2D) XY model and find the topological current relating to the order parameter field. 
A scalar field, -^j, is introduced through the topological current theory. By solving the scalar field, 
the interaction energy of vortices in a 2D XY model is revisited. We study the dynamical evolution 
of vortices and present the branch conditions for generating, annihilating, crossing, splitting and 
merging of vortices. During the growth or annihilation of vortices, the dynamical scaling law of 
relevant length in a 2D XY model, £,{t) oc (f — t*)^^^ , is obtained in the neighborhood of the limit 
point, given the dynamic exponent z — 2. This dynamical scaling behavior is consistent with 
renormalization group theory, numerical simulations, and experimental results. Furthermore, it 
is found that during the crossing, splitting and merging of vortices, the dynamical scaling law of 
relevant length is ^(t) oc (t — t*). However, if vortices are at rest during splitting or merging, the 
dynamical scaling law of relevant length is a constat. 

PACS numbers: 47.32.C-, 61.72.Cc, 64.70.qj 



I. INTRODUCTION 

Vortices play an important role in understanding a va- 
riety of problems in physics. In the 1970s, Kosterlitz and 
Thouless constructed a detailed and complete theory of 
2D systems [l| . The KT phase transition theory predicts 
that vortices pair unbinding will lead to a second-order 
transition in 2D systems, such as superfluid films, super- 
conductors, and XY models 0, Q. In a 2D XY model, 
there exist meta-stable states corresponding to vortices, 
which are closely bound in pairs below a critical tempera- 
ture. Above the critical temperature, the paired vortices 
unbind and become free. Vortices in a 2D XY model dis- 
rupt the spin alignments even at large distances and cor- 
respond to singularities of the order-parameter field Q . 

In 2D XY model, vortices are topologically stable con- 
figurations. It is found that the high temperature dis- 
order phase with an exponential correlation is a result 
of the formation of vortices. The critical temperature at 
which the KT transition occurs is, in fact, that at which 
vortex generation becomes thermodynamically favorable. 
At temperatures below this, the system has a power law 
correlation. In a low temperature phase, vortices appear 
in a small density of tightly bound dipole pairs. With 
an increase in temperature, vortex pairs dissociate and 
become free in the disorder phase. The evolution of vor- 
tices plays an important role in the KT transition of the 
2D XY model. 

There has been progress on the study of the defects 
associated with an n-component vector order parameter 
field, (/)(r, t). For the scalar case (n=l), the defects are 
domain walls, which are points of the spatial dimension- 
ality d=l, lines for d=2, planes for d=3, etc. More gener- 



* Email: ychenSlzu. edu. en 



ally, for n=d, one has point defects. This leads to vortices 
for n=d=2. It is interesting to consider the appropri- 
ate form for the point defect densities when expressed in 
terms of the vector order parameter field. This has been 
carried out by Halperin [5|, and exploited by Liu and 
Mazenko however, their analyses are incomplete 
In a 2D system, a gauge field-theoretical formalism has 
been developed by Kleinert Furthermore, the gauge 
theory of topological quantum melting in a 2 -|- 1 dimen- 
sional Bose system was developed by Nussinov et al, and 
the superfluidity and superconductivity can arise in a 
strict quantum field-theoretical setting [9]. 

A topological field theory for topological defects has 
been developed by Duan et al [l^. By using the 0- 
mapping method and topological current theory, the evo- 
lution of the topological defect that relates to singular- 
ities of the order-parameter field, such as the vortex in 
BEC [ll| and superconductivity was studied. In 
this paper, we will discuss the topological quantization 
and evolution of vortex in two-dimensional XY model. 
We introduce a scalar field, ip, through the topological 
current theory. By solving the scalar field, the Hamilto- 
nian of the interaction of vortices in the 2D XY model is 
revisited. 

More recently, the dynamical behavior of the 2D XY 
model following a quench of the temperature to below 
the KT critical temperature has been studied [l3|. Dur- 
ing the quench of a dynamic 2D XY model, aging phe- 
nomena and dynamic scaling behavior become extremely 
important [l^ . The assumption of dynamical scaling for 
the predicted asymptotic growth law the characteristic 
length, ^(t) (X {t/lnty/^ jl5j, which is also character- 
istic of the spacing between defects. However, scaling 
violations were reported. According to an expansion in 
t = A — d using standard field-theoretic renormalization 
group theory, it can be shown that the growing length 
is given as £,{t) cx t^^^ for large t, where z is the criti- 
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cal exponent for equilibrium critical dynamics pi| . This 
standard theoretical approach also shows that the result 
^(t) oc t^/^ is independent of the initial conditions. How- 
ever, this renormalization group method does not involve 
the effects of topological defects, such as vortices in the 
2D XY model. In the recent work, it was shown that, for 
the specific case of the 2D XY model, £^{t) oc {t/lnty/^ if 
free vortices are present, while S,{t) oc t^^^ if there are no 
free vortices present in the initial state • 

Through our topological current theory for the 2D 
XY model, the dynamic of vortices is studied, and the 
branch conditions for generating, annihilating, crossing, 
splitting and merging of vortices are given. During the 
growth or annihilation of vortices, the dynamical scal- 
ing law of relevant length in the 2D XY model, which 
is ^(i) (X {t — t*y^^, can be obtained in the neighbor- 
hood of the limit point. This indicates that vortices in 
the 2D XY model are the source of the scaling viola- 
tions. This dynamical scaling behavior is consistent with 
numerical simulations [3] and experimental results [l^ . 
Furthermore, we have also found that during the crossing, 
splitting, and merging of vortices, the dynamical scal- 
ing law of relevant length in the 2D XY model, which is 
£,{t) oc {t — t*). It has been shown that why the rel- 
evant length of vortices scales with time as t^^/^ for 
small values of times and then deviate toward a linear 
dependence in nematic liquid-crystal experimental ob- 
servations [23|. However, if vortices are at rest during 
splitting and merging, the dynamical scaling law of rel- 
evant length is f = const. Moreover, it is worthwhile to 
note that the dynamical scaling law of vortices, which 
was deduced from the topological current theory, only 
depends on topological properties of the order parameter 
field. 

The organization of this paper is as follows. In Sec. 
II, we describe the 2D XY model and the (/)-mapping 
method, and topological current theory is discussed. The 
Hamiltonian of vortices in the 2D XY model is revisited 
in section HI. The evolution of vortices and dynamical 
scaling discussed in section IV. Finally, in section V, we 
summarize our results. 



II. TOPOLOGICAL CURRENT IN THE 2D XY 
MODEL 

The 2D XY model is a system of spins confined to 
rotate in the plane of the lattice. The Hamiltonian of 
the system is give as: 



(1) 



<ij> 



where J is the strength of the nearest-neighbor interac- 
tion. 6i denotes the angle of the spin on site i with re- 
spect to arbitrary polar direction in the 2D vector space 
containing spins. 

By using the continuum limit, we can approximate 



the Taylor expansion. We can express partial derivatives 
through 9i — 6j = 89 for the two sites i and j . This leads 
to the continuum Hamiltonian: 

H = Ho-Eo = ^ J dr{Vef, (2) 

where Eq = 2JN is the energy of the completely aligned 
ground state of N spins. 

Kosterlitz and Thouless suggested that the disorder- 
ing is caused by topological defects, such as vortices in 
two-dimensional XY model, which are characterized by 
a mapping from some loop F in real space onto the order 
parameter space. For two-dimensional XY model, this 
implies 



d0= (^Veds = 2T:W,W^0,±l,. 



,n, 



(3) 



where W is the winding number of the vortex. By in- 
troducing a 2D XY model, the local order parameter is 
defined as: 



= = (4) 

Quantized vortices are topological objects associated 
with topological properties of the order parameter 5*. It 
is worth noting that the phase of the order parameter 
is undefined at the vortex core. In other words, vortices 
correspond to singularities of the order-parameter field. 
We define the unit vector field n as 



= V^^, a= 1,2. 



(5) 



where n°ri° — 1. From the unit vector, n{x), we can 
construct a topological current of the order parameter 
field in the 2D XY model, which carries the topological 
information of (/>(a;): 

jk ^ _L ^^^.^ A; = 0, 1, 2. (6) 

ZTT 

We will see that this current does not vanish at the zero 
point of (t>{x) or the singularities of the unit field n(x). 
By using the 2D topological current theorem, Eq. ^ can 
be rewritten in the compact form; 



J 



(7) 



It is the important relation between the (5— like topologi- 
cal current of vortices and the order parameter, ^P, where 
J^{(t>/x) is the vector Jacobian of </): 

J" = \e^'''eab^^r^Jcj,\ (8) 

or 



x 



by the first two terms 1 
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According to the implicit function theorem, the Jaco- 
bian's determinant can be given as: 



(9) 



The sohitions of the zero point of 4>{x) can be generally 
expressed as: 



X^Xi{t), y = yi{t), Z = l,2,. 



(10) 



which represent N zero points, zi{t) (1=1,2,. . . ,N), or a 
world line of N vertices in space-time. 

With the 5-function theory, (5^ (</>), can be expanded as: 



N 



1=1 



\JWx) 



(11) 



where the positive integer, /3;, is called the Hopf index of 
map a; — > 0. The meaning of /3i is that when the point 
r covers the neighborhood of the zero, zl, once the vec- 
tor field, (j), covers the corresponding region for /3/ times. 
Using the implicit function theorem and the definition of 
the vector Jacobian (Eq. [5]), we can find the velocity of 
the l-th defect, 



vi 



dzi 



J[<i>lx) 



Jicj)/x) 

j(^) = [Ji(0/x),j2(0/x)] (12) 

The spatial and temporal components of the defect 
current, j'', can be written as the form of the current 
and the density of a system of N classical point particles 
moving in a (2 -f l)-dimensional space-time, 



N 



Y,Pimvi5\r-m)) 



1=1 



N 
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1=1 



where rji is the Brouwer degree, 
J{c^/x) 



\JWx)\ 



±1. 



(13) 



(14) 



(15) 



It can clearly be seen that Eq. (I13|) shows the movement 
of vertices. The topological charge of vertices in the 2D 
XY model are conserved: 



dp 
dt 



= 0. 



(16) 



In addition, there is a constraint of charge neutrality: 



1 ^ 

■' 1=1 



(17) 



which indicates that vortices in the 2D XY model appear 
in pairs. 



III. THE HAMILTONIAN OF VORTICES IN 
THE 2D XY MODEL REVISITED 

In analogy to the velocity field in a superfluid, the dis- 
tortion field, u = V0, carries the topological information 
of vortices in the 2D XY model. By using Eq. ([4]) and 
(l5|), we can prove that 

and the vorticity is given as: 

V xu = B,{e''''eabdjn''dkn^) (18) 

where ek{k — 1, 2, 3) are the base vectors in the Cartesian 
coordinate system. Comparing Eq. p8)) and Eq. ([6]), we 
conclude that 



V X -u = 27rj°2 



(19) 



Therefore, in the 2D XY model, the vorticity of distortion 
field, V X u, can be expressed in terms of the topological 
current of the order parameter field. By comparing the 
5— like topological current Eq. ([7]) and Eq. (|19p . we have 
the important relation between vorticity and the order 
parameter field in the 2D XY model: 



' x 



(20) 



From Eq. ([20|) . we see that the vorticity of distortion 
field, V X M, does not vanish at the zero points of ^ . The 
location and direction of the ith vortex are determined 
by the ith singular point, i;(t), and the vector Jacobian, 
J{(j)/x), on zi{t), respectively. In the absence of vortic- 
ity, there are no zero values of the order parameter field, 
(5^(0) is zero, and Eq. (|20p becomes: 



V X Mo = 0, 



(21) 



which is the condition of irrotationality. Thus, Eq. (I^Dl) 
describes both the vortex-state and the irrotationality- 
state. To describe a collection of vortices at locations 
zi{t), we substitute Eq. (fM)) into Eq. (fT9|l. leading to: 



N 



V xu = 2ttJ2 Pirns'^ (r- zi{t))z. 



(22) 



1=1 



It can be see that Eq. represents N vortices that are 
charged with the topological charge, W. W is Pirn, which 
describes the inner topological structure of the vortex. 
To find the solution for Eq. (P^ . we can introduce a 
harmonious scalar field, ^ defining by: 



(23) 



In the absence of a vortex in the 2D XY model, we can 
see that the condition, didjO = djdi6, leads to V^V = 0. 
Thus, the 2D distortion can be written as, m = wq — V x 
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(zip), which contains two parts: the first part, uq = V(^, 
describe spin waves, while the second part is associated 
with vortices. By using 2D topological current theorem, 
it can be shown that: 



N 



(24) 



1=1 



The scalar field, ^p, behaves like the potential due to a 
set of charged particles. The solution of Eq. is given 
as: 



= ^/?^77/'"(l^-^l(^)l)• 



(25) 



1=1 



The continuum Hamiltonian of the vortices in Eq. ([2]) can 
be rewritten as: 

H-Y,Ec^^^- j{yi;f<f^=-^ yVVV,rf'x (26) 

where Ed is the core energy of the vortex. Substituting 
Eq. (p4|) and ([25|l into the Hamiltonian ([26|l . one can 
obtain: 



which is identical to the Hamiltonian of a 2D Coulomb 
gas with point charges of charge Wi = Piirji. The con- 
straint given by Eq. (|17[) is thus the constraint of charge 
neutrality. 



IV. THE EVOLUTION OF VORTICES IN THE 
2D XY MODEL 




The zero point of the order parameter field, ^P, which 
is associated with the locations of the core of vortices, 
plays an important role in describing the evolution of vor- 
tices in the 2D XY model. If the Jacobian determinant, 
J^{(t)/x) ^ 0, we will have the isolated solution of the 
zeros of the order parameter field in 2 + 1 dimensional 
space-time. But, when J^{(j>/x) = 0, the above results 
will change in some way and will lead to the branch pro- 
cess of vortices. We denote one of the vectors Jacobians 
at the zero points as {t*,zi). According to the values 
of the vector Jacobian at the zero points of the order 
parameter, there are usually two kinds of branch points, 
namely, the limit points and bifurcation points [IJ . Each 
kind corresponds to different cases of branch processes. 



A. The generation and annihilation of vortices 

We explore what will happen to vortices at the limit 
point {t*,zi). The limit points are determined by 



X 



0, jM- 

X 



= 0, JM - 

X 



^0 



^0. 



(28) 
(29) 



FIG. 1: Generating and annihilating of vortices pairs, (a) The 
origin of two vortices, (b) Two vortices annihilate in collision 
at the limit point. 



Considering the condition given by Eq. ([28|) and mak- 
ing use of the implicit function theorem, the solution of 
the zero points of (t>{x) in the neighborhood of the point 
{t*,zi) is given as: 



t = t{x), y^y{x), 
where t* = t{zl). In this case, one can see that: 



dx 



J\cplx) 



JWx) 



oo, 



or 



dt 
dx 



0. 



(30) 



(31) 



(32) 



The Taylor expansion of t 
{t*, zi) is given as: 

2dx^ 



t-t* 



t{x) at the limit points 
{x-zD^ (33) 



t'.Zi 
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FIG. 2: Two vortices collide with different directions of mo- 
tion at the bifurcation point in (2-|-l)-dimensional space-time. 
Two world lines intersect with different directions at the bi- 
furcation point; i.e., two vortices encounter at the bifurcation 
point. 



which is a parabola in the x-t plane. From this equation, 
we can obtain two solutions xi(t) and X2{t), which give 
two branch solutions (world lines of the vortices). If 



>0, 



we have the branch solutions for t > t* which are related 
to the origin of a dipole pair. Otherwise, we have the 
branch solutions for t < t* , which related to the annihi- 
lation of a dipole pair (Fig.[T|). 

Since the topological current is identically conserved, 
the topological charges of these two generated or anni- 
hilated vortices must be opposite of one another at the 
limit points, such as: 



Pivi + f^2m = 0. 



(34) 



This indicates that the vortices always generate and an- 
nihilate in pairs. From Eq. (j33p . one also obtains that 
the velocity of the vortices is infinite when they are an- 
nihilating, which agrees with the result introduced by 
Bray |22l |. Furthermore, a new result can be obtained 
that shows that the velocity of the vortices is infinite 
when they are generating, which is gained only from the 
topology of the order parameter field. For a limit point, 
it is required that, {4> / x)\f 7^ 0. A bifurcation point, 
on the other hand, must satisfy a more complex condi- 
tion. This case will be discussed in the following. 



FIG. 3: Vortices with the same direction of motion, (a) Two 
world lines tangentially contact; i.e., two vortices tangentially 
encounter at the bifurcation point, (b) Two world lines merge 
into one world line; i.e., two vortices merge into one vortex 
at the bifurcation point, (c) One world line resolves into two 
world lines; i.e., one vortex splits into two vortices at the 
bifurcation point. 



B. Encountering, Splitting and Merging of vortices 

Now, let us turn to consider the case in which the re- 
strictions on the zero point [t* ,zi) are given by: 



0, 



0,1,2, 



(35) 



which imply an important fact that the function relation- 
ship between t and x or y is not unique in the neighbor- 
hood of the bifurcation point {t* ,zi). This fact is easily 
seen from 



dx _ J^i4>/x) 
Itt ~ J{(i)/x) 



dy^ J\cj,/x) 
dt J{(t>/x) 



(36) 



which, under Eq. (j35p . directly shows the indefiniteness 
of the direction of the integral curve of Eq. ([55)1 at (i* , ij). 
For this reason, the point (t*,i/) is called a bifurcation 
point of the orientation order parameter. 

As we know, at the bifurcation point, (t* , zl), the rank 
of the Jacobian matrix, [d(f)/dx], is 1 in the 2D vector 
order parameter. With the aim of finding the different 
directions of all branch curves at the bifurcation point, 
we assume: 



dy 



^0. 



(37) 
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From the implicit function theorem, there is one function 
relation: 



y = f{x,t). 



(38) 



According to the (/)-mapping theory, the Taylor expansion 
of the solution of the zeros of the order parameter field 
in the neighborhood of (t*, zi) can be expressed as: 

A{x-zlf + 2B{x-zl){t-t*) + C{t-t*f + - • • = 0, (39) 

which leads to 



and 



C\^) +2B^ + A = 0, 
.ax J ax 



(40) 



(41) 



where A, B, and C are constants determined by the order 
parameter. The solutions of Eq. (|40l) or Eq. ([4T|) give 
different directions for the branch curves (world line of 
the vortices) at the bifurcation point. There are four 
possible cases which demonstrate the physical meaning 
of the bifurcation points. 

Case 1 (A ^ 0): For A 4(5^ - AC) > 0, from 
Eq. (|40|) we get two different motion directions of the 
core of the vortex given as: 



dx 

It 



-B ± - AC 

A ' 



(42) 



which is shown in Fig. [21 where two world lines of two 
vortices intersect with different directions at the bifurca- 
tion point. This shows that two vortices encounter and 
then depart at the bifurcation point. 

Case 2 {A ^ 0): For A = 4(5^ - AC) = 0, form 
Eq. (j40p . we obtain only one motion direction of the core 
of the vortex given as: 



dx 



1,2 



B 
7' 



(43) 



which includes three important cases (Fig. [3]). These 
three cases are: (i) Two world lines tangentially contact; 
i.e., two vortices tangentially encounter at the bifurcation 
point, (ii) Two world lines merge into one world line; 
i.e., two vortices merge into one vortex at the bifurcation 
point, (iii) One world line resolves into two world lines; 
i.e., one vortex splits into two vortices at the bifurcation 
point. 

Case 3 (^ = 0, C 7^ 0): For A = A{B^ - AC) = 0, from 
Eq-dlOl) we have: 



dt 
dx 



-B±^/W~AC 
C 



0, 



2B 
'C' 



(44) 



There are two important cases (Fig. U]): (i) One world 
line resolves into three world lines; i.e., one vortex splits 





FIG. 4: Two important cases of Eq. (|44)) . (a) One world line 
resolves into three world lines; i.e., one vortex resolves into 
three vortices at the bifurcation point, (b) Three world lines 
merge into one world line; i.e., three vortices merge into one 
vortex at the bifurcation point. 



into three vortices at the bifurcation point, (ii) Three 
world line merge into one world line; i.e., three vortices 
merge into one vortex at the bifurcation point. 

Case 4 (A=C=0): Equations (gHl) and (gT]) give respec- 
tively: 



dx 



0, 



dt 
dx 



0. 



(45) 



This case shows that two world lines intersect normally at 
the bifurcation point, which is similar to case 3. It is no 
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surprise that both parts of Eq. ([45]) are correct because 
they give the slope coefficients of two different curves at 
the same point {t*,zi). 

The remaining components of dy/dt can be calculated 
from a function relation of Eq. 



dy 
dt 



df dx df 
dx dt dt 



(46) 



The above solutions reveal the evolution of the vor- 
tices. The topological structure of the vortices is detailed 
in the neighborhood of the bifurcation points of the or- 
der parameter field. Besides the encountering of vor- 
tices, i.e., two vortices encountered at and then depart 
from the bifurcation point along different branch curves 
(Fig. [2] and Fig. [3Ka)), splitting and merging of vortices 
are also included. When multicharged vortices pass the 
bifurcation point, it may split into several vortices along 
different branch curves (Fig. [31(c), Fig. Hj^a)). On the 
other hand, several vortices can merge into one vortex at 
the bifurcation point (Fig. [3Kb), Fig. [H^c)). As before, 
since the topological current of the vortices is identically 
conserved, the sum of topological charges of the final vor- 
tices must be equal to that of the initial vortices at the 
bifurcation point, which is given as (for fixed 1): 



f 



(47) 



This indicates that vortices with a higher value of Burg- 
ers vector can evolve to the lower value of Burgers vector, 
or that vortices with a lower value of Burgers vector can 
evolve to a the higher value of Burgers vector through the 
bifurcation process. Furthermore, we see that the gen- 
eration, annihilation, and bifurcation of vortices are not 
gradual changes, but that they start at a critical value 
of a parameter, i.e., a sudden change. It is important to 
note that further bifurcations are possible during the evo- 
lution of vortices besides the cases studied in this work. 
It is necessary need to assume that the terms in the Tay- 
lor series considered above vanish and that the expansion 
of the field is dominated by higher order terms. 



Dynamical scaling law of vortices in the 2D XY 
model 



In the neighborhood of the limit point, we denote the 
scale length / = Vx. The growth velocity or annihilation 
velocity of vortices is v = I /\/t [5^] . From Eq. [33)1 , one 
can obtain the scaling law as: 



V a (t-r)-i/2. 



(48) 



It can be seen that Ek cx {t - 1*)^^ . In the 2D XY 
model, the growth or annihilation is parameterized in 
terms of a relevant characteristic length ^ (i ) , which is also 
characteristic of the mean vortex-antivortex separation 



distance. From Eq. ([48]) , we find that the relevant length, 
£,{t), obeys the following: 



^{t)r^{t-t*y/^ 



(49) 



Equation is the dynamic scaling law of the vortex- 
antivortex pairs. The low temperature equilibrium phase 
has essentially no vortex pairs, and ^(t) is infinite. 
The relationship between the relevant length and the 
metastable vortex density below the critical temperature 



IS given as: 



(50) 



then the number of vortices that satisfies the power 
law N oc t~^. The dynamic scaling law of the vortex- 
antivortex pairs is consistent with renormalization group 
theory and results from recent numerical simulations us- 
ing Langevin equations and Monte Carlo methods [TtI . 
|l8|. The specially prepared nematic liquid-crystal sys- 
tem [l^, which is developed, exhibits 2D XY behav- 
ior and can be utilized to obtains the dynamical scal- 
ing behavior of vortices (Fig. [5]). This result is consis- 
tent with the expansion in e = 4 — d using standard 
field-theoretic renormalization group theory. This stan- 
dard theory shows that the growing length, $,{t) cx t^^^ 
for large i, where z is the critical exponent for equilib- 
rium critical dynamics [16]. This standard theoretical 
approach also shows that the result ^(t) cx t^^^ is inde- 
pendent of the initial conditions. However, this renor- 
malization group method does not involve the effects of 
topological defects, such as vortices in the 2D XY model. 

From the assumptions of dynamical scaling law, it has 
been predicted that the asymptotic growth law of the 
characteristic length is given by ^(<) cx {t/lnty^"^ p^ : 
however, scaling violations were reported. In a recent 
work, it was shown that the specific case of the 2D XY 
model is given as ^(t) cx (t/lntY^^ if free vortices are 
present, while ^(t) cx t^^^ if there are no free vortices 
present in the initial state [13]. Since the topological 
current is identically conserved, there are no free vortices 
during the growth or annihilation process. 

In the neighborhood of the bifurcation point, we denote 
scale length Vx — /. From Eqs. ([42l) - (|44|) . we can then 
obtain the growth or annihilation velocity of the vortices, 
which is given as 

V cx const. (51) 
The approximation asymptotic relation of is then 



e(t)cx (t-r). 



(52) 



It can be seen that this is the reason why the relevant 
length of the vortices scales with time as t~^^'^ for small 
values of t and then deviates towards a linear dependence 
in nematic liquid-crystal experimental observations [13]. 
From Eq. one can obtain 



^(t) = const, V = 0. 



(53) 



8 



1 000 


. ' " 1 1 ' ' ' 


1 




100 
















10 


: • Pargcllis's experiments 








: A Loft's simulation (T=5.0->0.9) 


□ ^\ 






■ □ Loft's simulation (T=5.0->0.8) 




1 


■ ■ ■ 






1 10 


100 






Time (t) 







FIG. 5: Number of defects versus time. Experimental and 
simulation results are shown. The results of simulation were 
the average of 20 quenches [l^ . for which T* = 5.0 initially 
and is suddenly reduced to T* = 0.8 (□) and T* = 0.9 (A). 
The experimental results, which are obtained by Pargellis et 
al [3, were averaged over nine experimental runs (•). Solid 
lines are the expected A'' cx t^^ scaling. 



It can be seen that that vortices are relatively at rest 
when ^(t) — const. 

During the growth or annihilation of vortices, the dy- 
namical scaling law of the relevant length in the 2D 
XY model is given as ^{t) (x (t ~ t*y^^. This indi- 
cates that the vortices in the 2D XY model are the 
source of the scaling violations. During the crossing, 
splitting and merging of vortices, the dynamical scaling 
law of relevant length in the 2D XY model is given as 
^{t) (X {t — t*). However, if vortices are at rest during 
merging or splitting, the dynamical scaling law of rele- 
vant length is ^{t) = const. Moreover, it is worth noting 
that the dynamical scaling law of vortices depends only 
on topological properties of the order parameter field. 



V. CONCLUSIONS 

In summary, we have studied the inner structure and 
evolution of vortices in the 2D XY model by making using 
of the 0-mapping topological current theory. A scalar 
field, tp, is introduced through the topological current 
theory. By solving the scalar field, the interaction energy 
of vortices in the 2D XY model is revisited. 

By using the 0-mapping current theory, the densities 
of vortices in terms of the order parameter field in the 2D 
XY model are obtained directly from the definition of the 
topological charges of the vortices. The inner topological 
structure of the charge of the vortices, which is charac- 
terized by the Hopf index and the Brouwer degree, was 
obtained. By using the 2D topological current theorem, 
the Hamiltonian of the interaction vortices in the 2D XY 



model was revisited. This result is identical to the Hamil- 
tonian of a 2D Coulomb gas with point charges of charge 
= At?,. 

Furthermore, we have studied the evolution of vortices 
in the 2D XY model and concluded that there are cru- 
cial cases of branch processes in the evolution of vor- 
tices when J^{(j)/x) = and (f) = 0, i.e., ry; is indefinite. 
It is important to note that according to the Landau- 
Ginzburg approach, there is instability in the 2D XY 
model at the zero point of the order parameter fields in 
the low temperature phase. This indicates that vortices 
are unstable in the low temperature phase. In the high 
temperature phase, the zero point of the order parameter 
field is stable. Due to the branch condition, the vortices 
generate or annihilate at the limit points and encounter, 
split, or merge at the bifurcation points of the order pa- 
rameter field. This result also indicates that the velocity 
of the vortices is infinite when they are being annihilated 
or generated, which is obtained only from the topologi- 
cal properties of the order parameter field in the 2D XY 
model. The scaling law of relevant length, which is given 
as R{t) oc (t — t*y/^, can be obtained in the neighbor- 
hood of the bifurcation point during the growth or an- 
nihilation of vortices. During the crossing, splitting and 
merging of vortices, the dynamical scaling law of relevant 
length is ^{t) cx (t — t*). However, if vortices are at rest 
during merging or splitting, the dynamical scaling law of 
relevant length is f = const. The dynamical scaling law 
of vortices only depends on topological properties of the 
order parameter field. 

APPENDIX A: TOPOLOGICAL CURRENT 
THEORY 

We study a two-component vector order parameter, 

(j) = (</>^,0^), over the base manifold, M (in this paper 
M = i?2(g)i?): 

^^| ^ (j)^ + i4? , c/}" = (j)''{x,y,t), a =1,2. (Al) 

Quantized defects are topological objects associated with 
topological properties of the order parameter 5*. It is 
worth noting that the phase of the order parameter is 
undefined at the vortex core. In other words, vortices 
correspond to singularities of the order-parameter field. 
Let us define the unit vector field, n, as: 

where, n°'n'^ — 1. From the unit vector, n{x), we can 
construct a topological current of the order parameter 
field in the 2D XY model, which carries the topological 
information of (j){x): 

= ^e''''eabd^n'^djn\ fc = 0, 1, 2. (A3) 

It can be seen that this current does not vanish only at 
the zero point of (j){x), or the singularities of the unit 
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field n(x). Obviously, the current given by Eq. (|A3p is 
identically conserved: 



dkf = 0. 



(A4) 



Suppose there is a vortex located at Zi . The topological 
charge of the defect is defined by the Gauss map, n : 



W{ct>,zi) 



1 

2^ 



eatn^dn^ 



(A5) 



Using Stokes' theorem in the exterior differential form, 
one can deduce that: 



1 

2^ 



(A6) 



If we define the Jacobian as: 



(B4) 



and, by virtue of the Laplacian relation in (j> space [24 



A^(ln(||<^||)) = 27r<52(^), 



where 



A, 



d d 



(B5) 



APPENDIX B: TWO-DIMENSIONAL 
TOPOLOGICAL CURRENT THEOREM 

By using the 0-mapping method, we can prove that the 
topological current given by Eq. (|A3p is a 5-like current: 



(Bl) 



From the above formula, it can be seen that the topolog- 
ical current does not vanish only at the zero points of the 
vector field (t>{x). 

Substituting Eq. \A2\ into Eq. (|A3p and considering 
that: 



dun 



^^HM 



(B2) 



we have 



/■ = i.,^.'=e.,9.<^"5} A_^(in(||0||)), (B3) 



is the two-dimensional Laplacian operator in cj) space, we 
obtain a (5-function like current. This current is given by: 



X 



(B6) 



From the above formula, one can see that the topological 
current does not vanish only at the zero points of the 
vector field (t>{x). Therefore, it is essential to investigate 
the solutions of <j){x) 0. 
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